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We study the optical, Raman, and ac Hall response of the doped Mott insulator within the dynamical 
mean-field theory (d = oo) for strongly correlated electron systems. The occurence of the isosbectic 
point in the optical conductivity is shown to be associated with the frequency dependence of the 
generalized charge susceptibility. We compute the Raman response, which probes the fluctuations 
of the "stress tensor", and show that the scattering is characterized by appreciable incoherent 
contributions. The calculated ac Hall constant and Hall angle also exhibit the isosbectic points. 
These results are also compared with those obtained for a non-FL metal in d — oo. The role of 
low-energy coherence (FL) or incoherence (non-FL) in determining the finite frequency response of 
strongly correlated metals in d — oo is discussed in detail. 



PACS numbers: 75.30.Mb, 74.80.-g, 71.55.Jv 
I. INTRODUCTION 

The celebrated Landau theory of the Fermi liquid 
(FL) [0 , has been the mainstay of the conventional the- 
ory of metals for over four decades. It has proved to be re- 
markably stable even when local correlations are strong, 
as has been observed in a variety of heavy fermion met- 
als. Recent theoretical work concentrating on the d = oo 
limit of lattice fermionic models has shown up the ro- 
bustness of the local Fermi liquid picture, even for large 
values of the interaction . 

The d = oo works have also clarified the conditions un- 
der which the metallic phase for a given lattice model is 
described by FL theory Q ; in the situation where symme- 
try breaking is suppressed, the FL metal survives as long 
as there is no ground state degeneracy, as in the Hub- 
bard model. In every case where the symmetry-unbroken 
ground state exhibits a degeneracy as a function of the 
model parameters, FL behavior is invalidated; examples 
arc the spinless Falicov-Kimball model Q and the two- 
channel Kondo model || in this limit. The crucial role 
of the transfer of spectral weight over wide energy scales 
across the insulator-metal transition has also been re- 
vealed; it has thereby become clear that the itinerant 
and the atomic aspects of the problem should be treated 
consistently on a common footing to obtain reliable an- 
swers. 

Optical conductivity is an insightful probe which can 
provide detailed information about the finite frequency 
charge dynamics of a correlated electronic fluid M. In 
recent years, much information about the non-Fermi liq- 
uid charge dynamics in the "normal" state of cuprate 
superconductors has been gleaned by careful optical mea- 
surements The inplane optical conductivity, a a b(uj) 
for hole doped cuprates shows a non-FL fall-off with u> 
at small u> << D (D is the bandwidth) ||, while the 
electron-doped cuprate N ' di_ x G 'e x CuO \- y shows behav- 
ior characteristic of a FL . A curious feature observed 
in these studies is the fact that the various <J a b{w) curves 



for different dopings cross at a single point (this is the 
isosbectic point in the literature) . It is interesting to ask 
for the underlying physics manifesting itself in such ob- 
servations. 

Hall measurements provide a detailed picture of fi- 
nite frequency charge excitations in an external magnetic 
field, and probe the nature of the transverse scattering 
processes in an electronic fluid . Measurements carried 
out on the cuprates show that the transverse scattering 
processes are characterized by anomalous temperature 
and frequency dependences, leading to the two- relaxation 
rate phenomenology j?j . It is well known that FL trans- 
port is generically characterized by a single scattering 
rate governing transport, and the appearance of two re- 
laxation rates is therefore cited as a striking manifesta- 
tion of the breakdown of FL ideas. 

Given the above, the theoretical problem of computing 
the ac conductivity tensor for a strongly correlated metal 
is interesting; this quantity encodes detailed information 
about the finite frequency charge excitations and their 
response to external electric and magnetic fields. The- 
oretically, the problem of computing the transport co- 
efficients for a strongly correlated fermionic system in a 
controlled way is a rather hard task. The problem is even 
harder in the case of magnetotransport; to compute the 
Hall conductivity tensor, one has to evaluate explicitly 
three-point functions. In finite spatial dimensions, vertex 
corrections, which may be important, cannot be evalu- 
ated satisfactorily in any controlled approximation. The 
above difficulties make it imperative to search for con- 
trolled approximation schemes where some of the above 
difficulties can be circumvented without sacrificing essen- 
tial correlation effects. 

The dynamical mean-field approximation (DMFA or 
d = oo) has proved to be a successful tool to investigate 
transport in strongly correlated systems in a controlled 
way ||. This is because the vertex corrections entering 
in the Bethe-Salpeter eqn for the two-particle propagator 
for the conductivity vanish rigorously in d — oo. To eval- 
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uate the conductivity tensor, one needs only to compute 
the fully interacting local self-energy of the given model, 
following which the Kubo formalism can be employed |j] . 
Given that DMFA captures the nontrivial local dynamics 
exactly, one expects that it provides an adequate physi- 
cal description in situations where local fluctuations are 
dominant. In fact, the dc resistivity and Hall effect, as 
well as the longitudinal ac conductivity for the Hubbard 
model, have already been considered in the literature, in 
the framework of the d = oo approximation [ fL0| . 

However, the harder problem of computing the fre- 
quency dependent Hall effect (Hall constant and angle) 
has not been studied to date in detail. Actually, Lange 0] 
has employed the Mori-Zwanzig projection formalism to 
study dc magnetotransport in the Hubbard model. The 
actual evaluation of the complicated equations is, how- 
ever, actually carried out in the Hubbard I approxima- 
tion. This is known to lead to spurious instabilities (like 
ferromagnetism, which is washed away when local quan- 
tum fluctuations are included). Moreover, only the dc 
Hall constant is evaluated explicitly. The Hubbard I ap- 
proximation does not correctly capture the transfer of 
high-energy spectral weight to low energy upon hole dop- 
ing, a feature characteristic of correlated systems, and 
so one expects that it will be inadequate when one at- 
tempts to look at the ac conductivity. Some of the de- 
ficiencies of the Hubbard I approximation maybe cured 
by the so-called Hubbard III approximation jO]; this is 
actually the exact solution of the d = oo Falicov-Kimball 
model However, the metallic phase is not a Fermi 

liquid, and actually describes the charge dynamics in a 
model with X-ray edge singularities at low energy [Q. 
More recently, Lange et al. |?J have computed the ac 
Hall constant using the IPT in d = oo. However, the 
ac Hall angle and the Raman response, which we com- 
pute, has not been considered in Ref. 0|. Additionally, 
we have also compared our results with those obtained 
for a non-FL metal JllJ in d = oo to clarify the role of 
low energy coherence (incoherence) in determining the ac 
response of correlated metals. The detailed doping and 
frequency dependence of the Hall constant and Hall an- 
gle for a strongly correlated FL in d = oo thus remains 
an open problem to date. 

In this paper, we attempt to fill in this gap by study- 
ing the ac Hall effect and its doping dependence in de- 
tail. We concentrate on the one-band Hubbard model; 
the simplest model exhibiting FL behavior in d = oo. 
We have recently studied the ac Hall response for the 
"simplified Hubbard model" exactly in this limit- as 
mentioned above, it describes the response of a non-FL 
metal. We make an in-depth study of all the properties 
which can be computed without making further approx- 
imations. Specifically, we study the optical conductivity 
tensor, the ac Hall constant and Hall angle, as well as the 
Raman intensity lineshape as a function of filling for the 
correlated FL metal. We also compare our results with 



those obtained for a non-FL metal in d = oo. 

II. MODEL AND COMPUTATIONAL DETAILS 

We start with the one-band Hubbard model, 

H = - f E (clc ja + h.c) + UJ2^mi (1) 

defined on a hypercubic lattice in d dimensions. In the 
d — > oo limit, this corresponds to a gaussian unperturbed 
DOS In this limit, the lattice model, eqn.(l), is 

mapped onto a one-channel Anderson impurity problem 
(SIAM) embedded selfconsistently in a dynamical "bath" 
function that encodes the dynamical information about 
the quantum nature of the problem. Solution of the HM 
(eqn.(l)) therefore requires a reliable way to solve the 
impurity problem Unfortunately, there is no exact 
analytical solution available for the SIAM, and one has 
to resort to schemes which give reliable answers, and 
agree well with "exact" methods, e.g, with exact diag- 
onalization In this paper, we use the iterated pertur- 
bation theory (IPT) away from half-filling |ll| to solve 
the SIAM. This technique uses the Friedel sum rule |l4| ] 
(equivalent to the Luttinger theorem) to ensure that the 
correct Fermi liquid behavior is recovered at low energy. 
By construction, it is also exact in the band and the 
atomic limits, and so is a reliable interpolation scheme 
that describes the full local dynamical spectrum for all 
U/t and band-fillings. Knowledge of the local dynamics 
in the SIAM enables us to compute the full local dynam- 
ics of the HM in d = oo. We do not repeat the features 
of this method here, but refer the interested reader to 
Ref. [H for details. 

Solution of the d = oo problem using the IPT yields 
the full local self energy S(w) and the Green function 

GM = l/JVE k G(k,w) = 1/iVEkN - *k - s(w)]- 1 
for the symmetry unbroken paramagnetic case, which we 
consider here. We have computed the local DOS for the 
HM, and have checked that all features of the strongly 
correlated Fermi liquid metal are reproduced in accor- 
dance with jl3j. The quadratic energy dependence of 
im£(ix>) at small energy, the collective Fermi liquid peak 
that shifts to lower energy with hole doping, as well as the 
transfer of spectral weight from the upper- to the (central 
FL resonance + lower Hubbard band) are all reproduced 
well. 



III. OPTICAL CONDUCTIVITY AND RAMAN 
LINESHAPE 

As mentioned above, knowledge of E(w) for the lattice 
problem is the only input required to calculate the optical 
conductivity tensor in d = oo. Since vertex corrections 
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drop out in the two-particle eqn. for the conductivity in 
this limit, the longitudinal optical conductivity is given 
by a simple bubble diagram involving the full interact- 
ing, local propagators. Using the well known eqn. for 

a xx (iw) = — [ p (e) y~] G{e, iv)G(e, iu + iv) (2) 

ILU J ^ 

we have computed a xx {u>) for a given choice U/D = 3.0 
and for different band- fillings 5 = (1 — n) (D is the ef- 
fective bandwidth of the non- interacting model). The 
results are shown in fig. |l]. 
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FIG. 1. The optical conductivity of the Hubbard model 
with U/D — 3.0 for various band-fillings; 8 — 0.1 (solid line), 
8 — 0.2 (dashed line), 8 — 0.3 (dot-dashed line). <j xx (uj) for 
the FKM (non-FL) for U/D = 3.0 and 8 = 0.1 (bold-dashed 
line) . 

We see that a xx (u>) shows three distinct features wor- 
thy of mention: (1) The low-energy part (lu < 0.25) 
is characteristic of a renormalized FL metal, and cor- 
relates well with the behavior of the central FL peak 
in the DOS. (2) However, ct xx (lu) begins to rise again 
around lu/D ~ 1, passes through a broad peak centered 
around U/2 (for n = I), and starts falling off at larger 
frequencies. The transfer of optical spectral weight from 
high- to low-energy upon hole doping is clearly exhib- 
ited, and is understandable in terms of the interplay be- 
tween the high energy (localized) part, which inhibits 
double occupancies, and the low energy (coherent) part 
of the spectrum, reflecting carrier itinerancy that in- 
creases with hole doping. (3) Most interestingly, the 
& xx {<+>) curves for various S all cross at a single point 
(lu ~ IAD); this is the isosbestic point cited in earlier 
work [Q . We have a partial understanding of this curious 
feature: Following Vollhardt Jig ], who investigated such 
features seen in the specific heat for the HM for different 



U , we ask for the reasons for the crossing of the different 
curves, and the width of the crossing region. Focussing 
on the first part of the question, it is easy to show that 
there must exist some to = lu c (u) for which the curves 
cross. First, consider the high-energy limit of our re- 
sults. In general, <j xx (lu) — — x"(lu)/lu. Now, as lu — » oo, 
E(w) = U 2 (n/2)(1 - n/2)/u, and G(uu) = 1/oj, so that 
<j xx (lu — * oo) = (n/2)(l — ri/2)/w, whereby da xx /dn > 
for lu — > oo for all n < 1. On the other hand, it is clear 
from fig. (I), as well as from the optical sum rule, that 
v xx {u) increases with increasing 5, i.e with decreasing n, 
at energies lu < IAD, so that we have da xx /dn < at 
small lu. This implies that cj xx (lu) curves for different S 
must cross at some lu c (u). This crossing occurs at a sin- 
gle point if Lu c (n) is independent of n. Our results show 
that this seems to be the case to a high accuracy, but 
we are unable to quantify this argument further. It is 
an interesting problem to inquire into the deeper reasons 
for such features in the HM (see Ref. [l5) for a thorough 
discussion of the crossing points in C P (U,T)). 
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FIG. 2. Raman intensity lineshape I xx (lu) for the Hubbard 
model (FL) with U / D = 3.0, and for various band-fillings; 
8 = 0.1 (solid line), 8 — 0.2 (dashed line), 8 = 0.3 (dot-dashed 
line). I xx (u)) for the FKM (non-FL) for U/D = 3.0 and 
8 = 0.1 (bold-dashed line). 

Knowledge of the longitudinal optical conductivity al- 
lows us to compute the Raman spectrum within the 
d = oo methodology. This is because, as long as o xx (lu) is 
completely determined by the renormalized bubble con- 
tribution (vertex corrections drop out, as in d = oo), the 
Raman intensity is directly related to the optical conduc- 
tivity [n6| via the eqn., 



I xx (uj) 



1 



—Rea xx (uj) 



(3) 



In fig. g, we show the Raman lineshapes calculated 
from the above eqn. for the same U /D and doping values. 
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Since all calculations are carried out at T = 0.01, I xx (oj) 
is finite at to = 0. A sharp peak at low energy, reminis- 
cent of the plasmon peak in continuum treatments Jl6| , 
is visible, followed by an incoherent response character- 
istic of Mott-Hubbard systems at higher energies. 
However, the low energy peak is more the characteristic 
of the coherent part of the particle-hole response at low 
energies in a strongly correlated Fermi liquid. Evidence 
in favor of this interpretation is provided by the fact that 
the sharp peak broadens out as the temperature is raised 
above Tjf, the lattice Kondo temperature, where it is 
destroyed by strong scattering off local moments (which 
exist for T > Tk, rendering the FL description invalid). 

It is interesting to compare the above results with those 
obtained for a non-Fermi liquid (NFL) metal in d = oo. 
In this limit, NFL behavior at T = is obtained when 
the local Kondo effect is suppressed; this is exactly what 
happens in the Falicov-Kimball model (FKM) in this 
limit 0. Since the |-spin hopping is zero in the FKM, 
the Kondo effect does not occur, and the metallic state 
is not a FL near half-filling (it is actually described by 
the local X-ray edge physics in d = oo). In the inset of 
figs. we show the real part of cr xx (u/) and the Raman 
intensity lineshape for the FKM with identical values of 
U/D = 3.0 and hole doping (S = 0.1). It is seen that 
a xx (ui) for the FKM falls off much more slowly (like 1/w) 
in comparison with the fast Drude-like fall off observed 
for the HM. As expected, the sharp low energy peak in 
I xx {w) is changed into a broad continuum, reflecting an 
incoherent low-energy response characteristic of a NFL 
metal. We also observe that I xx (io = 0) ^ at T = 
in the FKM case, as argued by Shastry et al. |l6fl . This 
corresponds to the fact that, in a non-FL metal (where 
the Green fn. has a branch cut rather than a pole struc- 
ture), the action of the kinetic energy, or the stress tensor 
(which has non- vanishing matrix elements between lower 
Hubbard band (LHB) states) does not create well-defined 
elementary excitations. This results in the broad contin- 
uumwith non- vanishing intensity in the FKM (NFL), in 
contrast to the sharp peak with vanishing intensity (at 
T = 0) in the HM (FL) observed above. Hence, the finite 
frequency response of correlated metals is determined by 
the presence of a pole (branch cut) structure in the single- 
particle Green fn. at low energy. Within the d = oo ideas 
used here, the low energy coherence is a manifestation of 
the collective lattice Kondo effect in the Hubbard model, 
while the incoherent response in the non-FL case is un- 
derstood in terms of the vanishing of the .[-spin hopping, 
and to the complete suppression of this Kondo scale. 



IV. AC HALL CONSTANT AND ANGLE 

As mentioned above, computation of the magneto- 
transport is an extremely delicate matter, since one has 
to evaluate three-point functions |Q to first order in the 



external magnetic vector potential A. Generally speak- 
ing, in a nearly-free electron picture, the Hall effect is 
determined by the vagaries (shape and size) of the Fermi 
surface. That such a correspondence cannot be made for 
strongly correlated metals was pointed out by Shastry et 
al. |l7j who showed that the Hall constant in a strongly 
correlated system is dominated by spectral weight far 
from the Fermi surface, and hence is independent of its 
shape. This suggests that the results should not be sensi- 
tive to the choice of the free DOS, allowing us to use the 
d = oo DOS for the hypercubic lattice. Calculations for 
a three-dimensional system, where the d = oo approach 
works quite well Q , can be carried out by replacing the 
d = oo DOS by a 3d DOS. 

To compute the Hall conductivity, we need to add a 
Peierls coupling term to the HM, where it enters via the 
hopping |17| . The Hamiltonian in a magnetic field is 



H 



<ij>,a 



Uj (A)(Cl C J(T + h.c) + U^ fHinu (4) 



where the hopping matrix elements are modified by a 
Peierls phase factor and are = exp(2in / (j) A • dl), 
where A is the vector potential and </>o = hc/e. The 
off-diagonal part of the conductivity involves the com- 
putation of a three-point function to first order in the 
external field, as mentioned before. Fortunately, a con- 
venient form has been worked out by Lange Q , so we use 
the approach developed there. Explicitly, after a some- 
what tedious calculation, the imaginary part of a xy (uj) is 
given by 



' xy 



^xy 

1 

Ul 



+ oo 



dep {e)e / duJiduj2A(e, u>i)A(e, u> 2 ) 
F(c,uji;uj) - F(e, uj 2 ; v) 



UJl - UJ 2 



+ (fjj — » —u>) 



(5) 



where 



F(e,w;wi) = A(e,wi - w)[/(wi) - - u)] (6) 

and A{e, lo) = —Im[uj — e — S(w)] _1 /7r is the s.p spectral 
function in d — oo. 

Given the s.p spectral fn. in d = oo, o~ xy "(ui) is com- 
puted from the above, and the corresponding real part 
is obtained from a Kramers-Kronig transform. The Hall 
constant and Hall angle are obtained directly as, 



and, 



cotdniuj) 



0~xy{w) 



(7) 



(8) 



We describe the results obtained for the HM with 
U/D = 3.0 and S = 0.1,0.2,0.3. In fig. g, we show the 
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real part of the ac Hall constant, R' H (cu), with the cor- 
responding imaginary part in the inset. Since we work 
at a small T = 0.01, the u> = part of R' H is finite and 
negative, in agreement with |hJ (R'^ is identically zero 
at u> = 0). A sharp low energy peak, whose strength 
grows with increasing doping, is clearly seen. It is asso- 
ciated with intraband transitions within the lower Hub- 
bard band. This is related to the fact that the kinetic 
energy has matrix elements between the lower Hubbard 
band (LHB) states, which correspond to quasicoherent 
processes with weight increasing with S. Interestingly 
enough, the isosbectic behavior shown in the optical con- 
ductivity is also shown up in R' h (lo). The corresponding 
results for the FKM (NFL) are shown with the bold- 
dashed line in fig. ||. Two features are worthy of mention. 
(1) R' H (u> = 0) > 0, in contrast to what is observed for 
the HM (FL), where it is negative, and (2) the low energy 
peak in the HM is replaced by a smoothly decreasing Hall 
response as expected. The underlying reason for this is 
again that in the FKM, the kinetic energy (stress tensor), 
which connects LHB states, does not create well-defined 
elementary excitations (since the Green fn. is totally in- 
coherent). 




FIG. 3. AC Hall constant of the Hubbard model with 
U/D = 3.0 for various band-fillings; 6 = 0.1 (solid line), 
8 — 0.2 (dashed line), 8 — 0.3 (dot-dashed line). Rh(u>) for 
the FKM (non-FL) for U/D = 3.0 and 8 = 0.1 (bold-dashed 
line) . 

In fig. ||, we show the real and imaginary parts of the ac 
Hall angle computed from the above eqn. It is clear that 
the strong frequency dependence of cotOniw) found near 
half-filling (6 = 0.1) is weakened on increasing hole dop- 
ing. The imaginary part is zero for lo = 0. This is again 
in marked contrast to the results of a similar calculation 
done for the NFL case (FK model), results for which are 
shown with the bold-dashed line in fig. [|. We have also 
computed the frequency dependent transverse scattering 



rate, T xy [ui) which goes like ui 2 in the FL case. Since the 
transport scattering rate as deduced from the Drude-like 
response of ct xx (uj) also varies quadratically in the d = oo 
Hubbard model, the longitudinal and transverse scatter- 
ing rates are described by a single relaxation rate, as they 
should be in a FL metal. This is in contrast to the re- 
sults obtained from a similar analysis for the FKM [jll| , 
where the longitudinal and transverse responses to an ap- 
plied Lorentz field are governed by qualitatively different 
timescales. 
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FIG. 4. AC Hall angle, cot0 H (uj) for the Hubbard model 
with U/D = 3.0 for various band-fillings; 8 — 0.1 (solid line), 
8 = 0.2 (dashed line), 8 — 0.3 (dot-dashed line). cotOniui) for 
the FKM (non-FL) for U/D = 3.0 and 8 = 0.1 (bold-dashed 
line) . 

An amusing feature of the present calculations done for 
the Hubbard model is the existence of isosbectic points 
in various calculated quantities, like a' xx (u>), Rh{u), and 
cotOni^), as a fn. of hole doping (see figs.(|l|-||)). It is an 
extremely interesting problem to examine this feature in 
more detail, perhaps in a manner analogous to Ref. [ p5| . 
We leave this issue for a more detailed investigation in 
future. 

A detailed investigation of the magneto-optical re- 
sponse for the d — oo Hubbard model, which in prin- 
ciple is computable from the conductivity tensor, will be 
reported separately. 



V. CONCLUSIONS 

In conclusion, we have determined the optical, Raman, 
and ac Hall response of the d — oo Hubbard model 
using the IPT off half-filling as a reliable approxima- 
tion iQ. We have identified interesting features in the 
results (isosbectic points mentioned above), and shown 
that the response of a strongly correlated Fermi liquid 
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metal is explainable in terms of the competition between 
the atomic and itinerant aspects inherent in the Hub- 
bard model. The evolution of the ac response with hole 
doping is controlled by the increasing weight of the qua- 
sicoherent processes (corresponding to transitions involv- 
ing only the lower Hubbard band states) relative to that 
of the high-energy incoherent features, whose weight di- 
minishes with progressive hole doping. Lastly, we have 
compared our results with those computed earlier by us 
for a non-FL metal in d = oo, and have discussed the role 
of low-energy coherence (or incoherence) in determining 
finite frequency response of strongly correlated metals. 
Given the success J2| of the d = oo methodology in un- 
derstanding aspects of the physics of three-dimensional 
transition-metal oxides, we believe that the calculation 
presented here in combination with the actual 3d band- 
structure (which is equivalent to using the LDA DOS, for 
e.g) can be more relevant for studying ac response of 3d 
TM oxides. 
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